Inelastic scattering and heating in a molecular spin pump 
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We consider a model for a spin field-effect molecular transistor, where a directed pure spin current 
is controlled by an external electric field. Inelastic scattering effects of such molecular device are 
discussed within a framework of full counting statistics for a multi-level molecular system. We pro- 
pose that the heating of the molecular junction can be controlled by external electric and magnetic 
fields. Characteristic features of the model are demonstrated by numerical calculations. 

PACS numbers: 85.65.+h 85.75.Hh 71.38.-k 73.20.Hb 



I. INTRODUCTION 

Fast development of experimental techniques allowing 
for miniaturization of electronic devices led to renewed 
interest in theoretical research in the area of quantum 
transport. In particular in molecular electronics commu- 
nity focus of the research was shifting from elastic (Lan- 
dauer) to inelastic charge transport through molecular 
junction, to its noise characteristics and heating. Spin 
flux as an alternative to charge current in electronic de- 
vices, where magnetic field provides an additional (to 
bias) mean of control, is studied by spintronicsJ^ Quan- 
tum ratchets^ and electric potential^"— were proposed 
as additional controls of the spin flux. Schemes for op- 
tical control of a spin trapped in quantum dots were re- 
ported in the literature^"— Recently a combination of 
spintronics with molecular electronics started to reveal it- 
self as molecular spintronics J^"— Theoretical schemes for 
spin pumps were considered in Refs. [l8l - l20l . and inelas- 
tic effects of spin transport through the junctions were 
reported for spin valves^i and for tunneling through a 
junction with embedded spinj^^i Shot noise of spin cur- 
rent was considered in Ref. [23l . 

Here we present a model for a molecular junction, con- 
sisting of a molecule between two normal metallic con- 
tacts, where pure spin current is controlled by an external 
electric field. The spin field-effect transistor is a gener- 
alization of a spin pump introduced in Ref. ITsl . Directed 
spin current seems to be a more convenient element for 
molecular spintronics device. As in the case of the spin 
pump, only pure spin current participates in the trans- 
port. We consider inelastic effects of the spin current re- 
vealed in the transport properties of the device. This con- 
sideration takes place as a part of a general approach of 
full counting statistics for multi-level molecular systems. 
In particular, we discuss current and zero frequency noise 
of the junction. 

Flexibility of molecules, and as a result well- 
pronoimced inelastic features in transport properties of 
molecular junctions, makes heating of a device an impor- 
tant issue in molecular electronics. Spin field-effect tran- 
sistors (FET) inherit this problem of the usual FETs, 
since spin current is also caused by electrons cross- 



ing molecule-contact interface which results in heating 
molecular device. Within the model we discuss spin cur- 
rent through molecular junction, and propose external 
electric and magnetic fields as possible controls capable 
of tuning molecular structure to diminish heating of the 
device. 

Section [ll] presents a model and outlines the method. 
Section [Till is devoted to full counting statistics of multi- 
level molecular systems. We describe general approach 
and introduce junction characteristics ~ spin and charge 
current and noise. Section IIVI is devoted to heating in 
spin field-effect transistors. In section |V] we present re- 
sults of numerical simulations. Section [VII concludes. 



II. MODEL 



We consider extension of the spin pump proposed in 
Ref. llSl to the junction situation. The spin pump of 
Ref. [l^ is a model of a two-level system coupled to one 
electrode. Application of magnetic field generates spin 
current in the contact due to spin-flip process taking 
place at the molecule. We propose a generalization where 
the spin flux is generated across the junction (spin cur- 
rent between two contacts). Note that the spin flux in 
our model is optically controlled, i.e. external electric 
rather than magnetic field controls spin flux. The junc- 
tion is composed of two molecules (represented by single 
levels 1 and 2 coupled to molecular vibrations uJi and 
IJJ2) attached to left L and right R normal metal elec- 
trodes, respectively. Each molecule is subject to a dc 
and ac pair of magnetic fields S,-''^'' and b\°''^'' [i — 1,2). 
Note that level dependent g-factors in nanowire quantum 
dots were recently reported in the literature;^ something 
that, here, is represented by the site dependent field. The 
molecules are weakly coupled by external source of light 
of particular frequency oje- The light can be used as a 
handle to switch on and off spin transport through the 
junction. The Hamiltonian of the system is 



H ^Ho + V, 



(1) 
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where 



-^0 = ^ Ekfika + ^ ( ^[£ia + M^{a^ + a\)\ 



keL,R;a 



i=l,2 ^ a 



J2 [VikdU^ka + H. 

^ {V2kdlcka + H.c}j, (2) 
keR-.a 



V ={VEe~'^^' + V^e'^^')J2id\j2. + H.c). (3) 



where {die) and (cka) are creation (annihilation) 
operators for corresponding state in the niolecule(s) and 
in the contact(s), n.^^ = dl^di^ and = cl^ha, a] (a.^) 
arc creation (annihilation) operator for vibration quanta 
of molecule i, a ^ ±1 is direction of spin projection, and 
where 



£iCT ~ Si — 7:9f^B^i 



We start by transforming the Hamiltonian ([T]) into the 
rotating frames of the magnetic fiel d^°i^^ 



SB{t) = - ^y^cr I flla + flko 

a \ keh 
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(5) 
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This transformation eliminates the time-dependence 
from the ac magnetic field terms, shifts the positions of 
molecular levels, and induces spin biases in the contacts, 
according to 



{B) 



£kc 



\ Efc+a^r /2 fcGi, 



(7) 
(8) 



Ek+auf'' /2 ke R 
The perturbation ((S]) takes the form 

(7 

In what follows we will put cjj = — cjj = to create 
a spin bias across the junction (see Fig. [T]for a sketch). 

Similarly, we put q^bB^''^^ = — g/isBj'^'^'' = wg. The 
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FIG. 1: (Color online) Sketch of an optically controlled spin 
field-effect transistor. 



interaction with the external optical field is taken into 
account within perturbation theory. Expansion of the 
evolution operator on the Keldysh contour up to second 
order in ([9]) in the rotating-wave approximation (RWA) 
leads to electronic self-energy due to interaction with ex- 
ternal optical field in the form (see Appendix [X] for de- 
tails) 



(4) = <^.,a'2|T4;pcosu;£;(t - t') 



(10) 



G'2^,2^(r,T')e-'""«(*-*'' 

Gi,,i,(r,r')e'""^(*-*') 



As a result of the transformation to the rotating frames 
and the RWA we obtain a time-independent (steady- 
state) description. 

The coupling between electronic and vibrational de- 
grees of freedom is treated below within two differ- 
ent approaches. The transport through the junction 
is considered in a basis obtained by a small polaron 
transformation^^ 



Sy 



Mr 



ai)nia 



(11) 
(12) 



i=l,2:a 



This transformation decouples the electronic and vibra- 
tional degrees of freedom on the molecule and dresses the 
molecular fermion operators 



difT^i 



exp 



K - ar) 



(13) 



We take the vibrational shift operators X into account 
within the usual Born-Oppenheimer-likc approximation, 
which allows us to introduce the Franck-Condon factors21 
multiplying the electronic GF 



G 



(t,t') EE ~i{Tcl4T)X,{T)dj„>{T')Xj{T')) 



^^{TA^r)d]^,{r')){T,X.,{T)X]{T')) 



In the description of the junction heating, we utilize 
the Born approximation^^ instead, which allows us to 
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keep the consideration simple. As usual we implement 
a non-crossing approximation, i.e. diagrams for electron 
transfer and interaction with vibrations do not cross (the 
processes do not happen simultaneously)!^ In this case 
Green function is obtained from the Dyson equation 

G..,,.'(r,r') = Gl°)^..,(r,r') (15) 

nil ,m2 CTi ,0-2 ^ ^ 
X ^mili ,m2CT2 (■^l ' ■^a) Gm2a2,ja' (^2 , t') 

where G^*'-' is Green function in the absence of electron- 
phonon interaction (Mi_2 = 0), and E^p'*) is self-energy 
due to electron-phonon interaction 

^'^l' (r, r') = z<5.,, |M,p A(r, r') G..,.., (r, r') (16) 

Here Di[T,T') = —i{Tca,i{T)a\[T')) is phonon Green func- 
tion. We treat it within a quasi-particle approximation.— 



III. FULL COUNTING STATISTICS 

The theoretical concept of full counting statistics was 
originally proposed by Levitov and Lesovik .'^^i'^^ The ap- 
proach was applied to the non-equilibrium Anderson im- 
purity model in Rcf. [33l Measurements of shot noise 
in molecular junctioi*^ prove possibility of experimental 
observation of moments beyond average current. This 
together with the flexibility of molecules, i.e. importance 
of inelastic effects in transport through molecular junc- 
tions, recently caused several theoretical investigations 
devoted to study of resonant level coupled to single vi- 
bration modelj^"— 

Here we discuss a generalization of the result derived 
by Gogolin and Komnik^ to a multilevel situation, and 
apply the expression to calculate the first (current) and 
second (noise) cumulants of the distribution. Follow- 
ing the derivation of Ref. [3^ for the case of multilevel 
molecule we obtain an expression for the derivative of 
the adiabatic potential in the form^^ 



t.) - -\ f Tr {E< (i.)e- G> - G>E< (i.)e-} . -\ f 
Tr<{S<(£;)e*^- [(i? - Hq - S— (i?)) [£>(i?)] (S - Hq + S++(i?)) + S<(i?) 
(i? - Ho + S++(£;)) [11<{E)] {E-Uo- ^-(E)) + S>(i?) 



(17) 



Here, i^tc) ^ counting field for the interface be- 
tween the molecule and contact K on the forward (back- 
ward) branch of the Keldysh contour, Xk = {X]^ — X~^)/2. 

' ', I-', '<', and '>' are causal, anti-causal, lesser, 

and greater projections. The trace runs over molecu- 
lar degrees of freedom. The electronic self-energy due 
to the coupling to contact K is denoted by Hk, and 
Sa is the total electronic self-energy dressed with the 
counting field A. In particular, within the non-crossing 



approximation's, 



S>-<(i?) 



E^^-^(^)^ 



''int.X 



(E), (18) 



where S^^^ ^ is the electronic self-energy due to interac- 
tions dressed with the counting field A. 

It seems difficult to obtain an expression for the log- 
arithm of the generating function (integral of pT]) over 
X]^) in the multi-level case. However, Eq. ([T7)) itself 
can be used to calculate cumulants. The time-averaged 
charge cumulant of order n due to charge transport 



through interfaces {Kj} (j = {1,2,..., n}) is 

(^"g)if„,...,A-i 



T 



(19) 



d 



d 



Here and below e ^ h ~ I. 

First cumulant yields the well-know expression for 
steady-state current^iiS 



Ik 



d 



d{iX 



K. 



(20) 



/ + 00 JTp 
^ —Tt[ii<{e)g>{e)-i:>{e)g<{e)\. 

Since no spin-fiip events are allowed on the metal- 
molecule interface, the electronic self-energy due to cou- 
pling to the contacts is diagonal in spin space, that is. 



Si^'(T,r')=5]|V^,.|V.(r,r'), 



(21) 



where i = liox K = L and i = 2 for = i?, whereas gka 
is the Green function (GF) of a free electron. As a result. 
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the charge current of electrons with spin a at interface 
K is given by 



e 
n 



-(-oo 



dE 
2^ 



(22) 



)<(i?) G> - 4^)>(i?) G< 

The spin and charge currents at interface K are 



As) 

'if 



(23) 
(24) 



respectively. 



The second cumulant yields an expression for the zero- 
frequency noise. In the non-interacting model, it reads^ 



Sk^kA^ = 0) = 



/ + 00 ip 
^ ^Tr{J;,,,K, [i:<AE)G>{E) + G<{E)i:>AE)]-iKAE)iKAE) 

(E) G++ (E) S>^ {E) G-{E) + (E) G" (iJ) S<^ {E) G++ (i?) 
S< ^ (i?) G> (E) S> ^ (i;) G< (i?) - S> ^ {E) G<{E)^<AE)G>{E)}, 

I 



(25) 



where ixiE) is a matrix of the energy- resolved current 
operator at interface K 

iK{E)^i:<{E)G>{E)-ll>{E)G<{E) (26) 

The spin and charge noise at interface K are expressed 

by 

^k'k ~ ^K't,K't — '^Sk^^kI + SK^,Ki, (27) 
^K^K ~ ^K'[,K^ + + SK^,Kii (28) 

where the spin resolved elements of the noise SKa.Ka' 
can be obtained from ([25]) by generalizing the counting 
field to include spin. 

Higher order cumulants for a multi-level molecule can 
be similarly derived from (|17|) in a, cumbersome but, 
straightforward manner. 



IV. HEATING 

Within the model (Fig. [1]), the charge current is com- 
pensated at each interface so that only spins arc trans- 
ferred through the junction. Nevertheless, since electrons 
cross the molecule-contact interfaces also in the spin field- 
effect transistor, the question of heating of the molecular 
vibrations is still important. To estimate the tempera- 
ture of the vibrations we extend an approximate scheme 
used by one of us in a previous publication4^ The essence 
of the approach is quasi-particle assumption used to de- 
scribe vibrational degrees of freedom (phonons) . In this 



case the only relevant characteristic of the vibration (be- 
sides frequency) is its average population, while actual 
nonequilibrium distribution may be disregarded (density 
of states is a delta- function). 

At steady-state infiux of energy through the molecule- 
contact interfaces, J™, is compensated by outflux, J""*, 
so that energy of the molecule does not change J = J"' — 
jout _ Q Since energy is carried by both electrons and 
phonons, we get 

Je + Jph = 0. (29) 

Within the quasi-particlc approximation the latter for a 
vibration w,; coupled to a thermal bath is^ 

Jph = ^^L0^ [NBEiuj,) - N,] , (30) 

where uji and iV^ (z = 1, 2) is a vibration and the corre- 
sponding average population of molecule i, respectively, 
whereas fi is the spectral function of the bath, and 
Nbe{^) = [e^" — is the Bose-Einstein distribution. 

The electronic energy flux for a system coupled to a 
set of baths B can be written as 

/ + 00 Jpi 
^ —ETr[j:<{E)G>{E)~j:>{E)G<{E)], 

(31) 

where is the electronic self-energy due to the coupling 
to bath B. In our model we want to estimate the temper- 
ature of two vibrations, thus, following the approximate 
procedure introduced in Rcf. |44j, we split the system into 
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two parts at the molecule-molecule interface. In this way 
we can speak about electronic energy flux in two parts of 
our system. Each part has two interfaces: one between 
the molecule and corresponding contact and the other 
between the molecules. The electronic self-energy due to 
the coupling to the contact is given by the standard ex- 
pression ()21[) . The self-energy on the molecule- molecule 
interface is given by pUj) . 

We obtain the electronic flux in the molecule i part 
of the system through the following procedure. We be- 
gin with the Hamiltonian H which is transformed to the 
rotating reference frames of the magnetic fields, Eq.®, 
however, without performing the small polaron transfor- 
mation. We, then, treat the electron-vibration coupling 
in the Born approximation^^. Within the non-crossing 
and quasi-particle approximations we accordingly obtain 
the electronic flux 

Je = u,Mf {[N, + - iV,4+') (32) 

where 

/ + 00 Jp 
^G<^^^{E)G>^^^iE±co.) (33) 

Using dSni) and ^ in 1^ yields 

The vibrational temperature is obtained from (|34p under 
assumption that the vibrations are populated according 
to the Bose-Einstein distribution. We note that a more 
physically motivated procedure to introduce the temper- 
ature of the molecular vibrations at non-equilibrium can 
be uscd)^ however, for demonstration purposes the sim- 
pler procedure described above suffices. 

Atomic cooling caused by sub-resonance optical exci- 
tation is a well-known effect4^ In molecular junction the 
analogous detuning may lead to cooling of the molecu- 
lar vibration. Here, we utilize the sub-resonance detun- 
ing of the electric and magnetic field frequencies from 
the energy difFercncc of the inter-molecular and contact- 
molecule electronic transitions as a mechanism to pump 
energy out of the molecular vibration. Naturally, cooling 
is most efficient when the detuning frequency coincides 
with the frequency of the vibration. Additional possibil- 
ities to cool the current carrying molecular junction are 
discussed in Ref. H^. 

V. NUMERICAL RESULTS 

Here we present numerical results in order to illustrate 
the transport properties of our model for the spin field- 
effect transistor sketched in Fig. [T] We use the level 
width r due to coupling to a contact as unit of energy. 
Unless explicitly specified otherwise, the parameters for 
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FIG. 2: (Color online) Elastic transport. Spin current vs. 
(a) ac magnetic field frequency lob (solid line, blue) and (b) 
electric field frequency ue- The latter shows results for Ve = 
(dotted line, black), 0.3 (dashed line, magenta), and 0.5 
(solid line, red) at ljb = 20. Also shown curve for Ve = 0.5 
at a;_B =0 (dash-dotted line, blue). Dashed line (black) in (a) 
shows charge current. See text for other parameters. 



calculations are temperature in the leads Te — 0.3, 
the positions of the levels in absence of external fields 
El = 62 = 0, escape rates to the contacts (wide-band 
approximation is employed) F^^^ = 1, strength of cou- 
pling to external de magnetic field uq = g^sB^'^'^^ ~ —10, 
strength of the coupling to the external ac magnetic field 
g/isB'"'^' = 0.01 and its frequency ujb = 20, strength of 
the coupling to the external electric field Ve = 0.5. The 
parameters for the molecular vibrations are indicated in 
each case separately. 

Figure[2]shows the spin current l''*^ in the elastic trans- 
port regime. In the calculations, the shift of the levels due 
to changes uylob is assumed to be compensated by the dc 
magnetic field, so that levels are set as in Fig. [T] and do 
not move. The spin current (solid line) vs. the external 
ac field frequency (spin bias) is plotted in Fig. [5^. Note 
that the electric field is an additional source of energy 
in the model, so that one will observe spin current even 
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FIG. 3: (Color online) Inelastic transport. Spin current 
vs. (a) ac magnetic field frequency ljb and (b) electric 
field fi'equency lue- The latter shows results for symmetric 
cji = a;2 = 3 (dashed line, red) and asymmetric uji = 3, 
[J2 = 4 (solid line, blue) cases. Electron-vibration coupling is 
Ml = M2 — 3. Other parameters are as in Fig. [S] 



without spin bias in the contacts. The sign of the spin 
current shows the direction of the spin-up flux (direction 
from left to right is chosen as positive). The flux switches 
around ws/F = 10, as is shown in Fig. [2^. The point of 
sign change is defined by a competition between the spin 
bias and the electric field mediated inter-molecular trans- 
fer. At =0 the lower spin levels are occupied, while 
the higher levels are empty (see Fig. [T]). In this regime 
the electric field facilitates the spin-up flux from right to 
left, which defines a negative sign of the spin flux. The 
spin bias grows with the frequency ljbi and the value 
ws/F = 10 the two processes cancel each other. When 
the spin bias grows further, the population of the lev- 
els changes. The higher levels become populated due 
to the increased bias, while the population of lower lev- 
els diminishes (the level goes above the corresponding 
spin-resolved chemical potential, see Fig. [1]). The posi- 
tions of the peaks are defined by the resonance condi- 
tion for inter-molecule electron transfer: the position of 
the molecular levels is kept fixed by the dc field adjust- 



ment, while the frequency of the transition changes with 
both electric field frequency loe and ac magnetic field fre- 
quency LdB, see Eq.([n]). The condition for the resonance 
iswBzbcjs = |£ict— £'2ct|- We will return to the question of 
the role the spin bias plays in our discussion of the noise 
and heating properties. Note also, that the charge cur- 
rent (dashed line in Fig. [2^) is identically zero (the flux 
of spin up electrons from left to right is compensated by 
the flux of spin down electrons from right to left). 

Fig. IDs shows the dependence of the spin current on 
frequency of the external electric field. Naturally, the 
dependence has a maximum at resonance. Shown are 
plots for magnetic field frequency ws in the first (dash- 
dotted line) and second (solid line) maximum of Fig. [2^. 
Also shown are curves for smaller (dashed line) and zero 
coupling to the electric field (dotted). In the absence of 
coupling no spin current is observed. This indicates that, 
for the parameter range chosen, we are not in the regime 
of the spin pump described in Ref. 

Figure [3] is inelastic analog of results presented in 
Fig. [21 Molecular vibrations are taken into account em- 
ploying small polaron transformation, and assuming sep- 
aration of timescales, so that introduction of the Franck- 
Condon factors becomes possible. Fig. [3^ shows spin cur- 
rent vs. frequency of ac magnetic field. Parameters for 
molecular vibrations are uji = lu2 = 2, electron-vibration 
coupling is Mi = M2 = 2. Other parameters are as 
in Fig. [2l An unusual form of the vibration sidebands 
is caused by resonance condition for intermolecular elec- 
tron transfer rather than resonance situation at molecule- 
contact interface. This makes vibrational characteristic 
different from the one presented in Ref. [2^ At the same 
time, vibrational structure observed in spin current vs. 
electric field frequency resembles such for charge current 
vs. bias plots. Figure shows two such characteristics: 
for symmetric uji = UJ2 = 3 (dashed line) and asymmetric 
wi < UJ2 {uji ~ 3 and uj2 = 4, solid line) cases. Presence 
of a higher frequency naturally leads to observation of 
less vibrational sidebands. 

Figure [4] displays the noise properties of the junction. 
We show the results for elastic transport only, leaving 
the study of inelastic noise properties for a future publi- 
cation. Fig. shows the spin (or charge) zero-frequency 
noise as function of the external electric field frequency. 
Note that while the charge current is identically zero, 
charge noise does exist. Shown are the results for two 
situations: the spin levels of each molecule are well re- 
solved (dashed line, red) vs. essential overlap between 
the two (solid line, blue). The first case corresponds 
to a situation with well pronounced resonance behavior. 
In such situations, the transmission probability at res- 
onance approaches unity, which leads to suppression of 
the noise4^ The latter case corresponds to a situation 
when the transmission probability is lower. In this case 
no noise suppression is observed at resonance. The effect 
is similar to the behavior of a molecular junction in sym- 
metric vs. asymmetric couplings at the two sides of the 
junction4^ 
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FIG. 4: (Color online) Noise properties, (a) Zero-frequency 
charge and spin noise for ri-^' = V (dashed line, red) and 
ri^' — 5 (solid line, blue), vs. electric field frequency uje- 
(b) Difference between spin and charge zero-frequency noise 
for giiBB^""^ = 0.01 (dash-dotted line, red), 0.05 (solid line, 
magenta), and 0.1 (dashed line, blue). See text for parame- 
ters. 



In contrast to charge noise, the spin noise strongly de- 
pends on spin-flip events within the system4^ Fig. |4)3 
shows the difference between the spin and charge zero- 
frequency noise at several magnitudes of the coupling to 
the external ac magnetic field. Note that at resonance, 
where the electron transport through the junction be- 
comes pronounced, the difference between the spin and 
charge components diminishes. 

As was mentioned already earlier, there is no net 
charge transport between the leads. Nevertheless, charge 
transfer between contacts and molecule does occur. 
Thus, the spin field-effect transistor may be heated. As 
was discussed in the literature (see e.g. Refs. I47ll50l ). 
charge transfer by itself does not necessarily lead to 
heating of the molecular device. Here, we demonstrate 
one such possibility. Cooling of a molecular vibration is 
caused by tuning the external field frequency out of res- 
onance (an analog of atomic optical cooling). Figure [5^ 




FIG. 5; (Color online) 'Temperature' of the first uji (dashed 
line, red) and second lu2 (solid line, blue) molecular vibra- 
tion vs. (a) ac magnetic field frequency ujb and (b) external 
electric field frequency uje- See text for parameters. 



shows the temperatures of the first wi (dashed line) and 
second ui2 (solid line) molecular vibrations as functions 
of the spin bias. The coupling to the dc magnetic field 
is ujQ = —4, whereas the frequency of the electric field 
is = 4. The molecular vibrations are taken to be 
cji = 1 and UJ2 = 0.5 with electron- vibration couplings 
Ml — M2 — 0.5. Other parameters of the calculation 
are introduced at the start of the section. The parame- 
ters (electric and magnetic field frequencies) are chosen 
in such way that the most probable (resonant) electronic 
tunneling through the junction has to go uphill in en- 
ergy at the contact-molecule and inter-molecule transi- 
tion steps. In this case, the tunneling energy difference 
required for resonant electron transition is taken from the 
molecular vibration, which leads to cooling of the device. 
Fig. [5)d shows similar behavior when the frequency of the 
electric field is used as a control instead. The molecular 
vibrations are wi = 0.5 and W2 = 0.2. Here cooling takes 
place at the molecule-molecule interface only and is not 
as effective. However, one can achieve a stronger signal 
in this case. The efficiency of the cooling at maximum 
spin current is an interesting question for future studies. 
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VI. CONCLUSION 

Spintronics is a quickly developing field of research 
(both experimental and theoretical). Here, we studied 
the spin current transport properties within a model for 
an all-electric controlled spin field-effect transistor. The 
model provides pure spin currents through a junction 
which consists of two spin-pumps (molecules under in- 
fluence of external dc magnetic fields)^ each attached to 
its normal metal contact. The two parts are coupled by 
an external electric field, which serves as a control for the 
spin flux. Within a full counting statistics approach to 
multi-level systems, we discuss clastic and inelastic spin 
flux and noise properties of the junction. External elec- 
tric and magnetic fields are indicated as possible controls 
of the spin current through the junction. The charge 
and spin zero-frequency noises are shown to be different 
when spin-flip processes within the junction become pro- 
nounced. Zero-frequency spin noise as function of the 
external electric field shows the same single to double 
peak structure transition when the tunneling probability 
approaches unity. Similar behavior was observed earlier 
for a model of charge field-effect transistor.— 

Problem of molecular junction heating recently dis- 
cussed in the literature in connection to charge transport, 
retains its importance also for spin molecular devices. 
Within the model, we discuss possibility to use external 
fields slightly de-tuned from the molecular resonances for 
cooling the molecular vibrations. This process is similar 
to optical cooling of atoms. We find that external ac 
magnetic fields may be effective cooling media. External 
electric fields also provides the effect. The efficiency of 
the cooling at maximum spin flux is an interesting ques- 
tion for future studies. 



Acknowledgments 

J.F. thanks UCSD for its hospitality during his visit 
in July 2009, and the Swedish Research Council and 
the Royal Swedish Academy of Sciences for flnancial 
support. M.G. gratefully acknowledges support by the 
UCSD (startup funds), the UC Academic Senate (re- 
search grant), and the U.S. -Israel Binational Science 
Foundation. 

Appendix A: Derivation of Eq. (|10p 

Here we derive an expression (|10p for the self-energy 
due to the coupling to the external electric fleld within 



the rotating-wave approximation. We start by partition- 
ing the total Hamiltonian of the system into a zero-order 

Hamiltonian Hq and perturbation V, Eq.®. The first 
is given by an analog of Eq.([2|) after transformation to 
rotating frames. It will be similar to ([2]) with molecular 
and contact states energies renormalized according to ^ 
and ([8]), respectively, and without the time dependence 
of magnetic field terms. For the moment, we disregard 
the coupling to molecular vibrations. Hamiltonian Hq 
defines the zero-order GFs of the system G*^"-*. We treat 
the interaction V by perturbation theory, i.e expanding 
the evolution operator on the Keldysh contour (interac- 
tion representation) exp[— i drV/ (r)] to second order in 
Ve- This leads to a Dyson type equation for the Green 
function 

G(r,r')=G(«)(r,r') (Al) 
+ jdTijdT2 G(")(r, n) S(^'(ti,T2) G(")(r2,T') 

with matrix elements for the self-energy given by 

4fJ„„..(n,r2) ^ SmAVEe-^--'- + y^e--*^) 
X (V^ije-'"^*^ + VEe"^"^^-) (A2) 

X G^°i^,^„jTi,r2)exp[i(-l)™a;B(o-iti - 0-2^2)] 

where m,n — 1,2 numerate molecule in the junction, 
whereas m means opposite of m. Application of the 
rotating wave approximation leaves only terms propor- 
tional to \Ve\^ and enforces cti = (T2, in (|A2[) . The re- 
sulting expression is presented in Eq. (jlO[) . 

In the case of the small polaron transformation, the 
coupling to molecular vibrations is treated by dressing 
the expression pO|) by Franck-Condon factors 

= '5™,„(5,,„,2|VE|2cOSW£(t - t') 

xG^.,^.(r,r')e'(-i''" (A3) 
X {T,XAT)Xi{r')){TcXAr')Xl{T)) 

When perturbation theory is employed, the GF in Eq. 
(|Aip is considered as a zero-order GF of the Born ap- 
proximation. 
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